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Abstract 

Rozowsky, Volkas and Wali pQ recently found interesting numerical 
solutions to the field equations for a gauged {7(1)®[/(1) scalar field model. 
Their solutions describe a reflection-symmetric domain wall with scalar 
fields and coupled gauge configurations that interpolate between constant 
magnetic fields on one side of the wall and exponentially decaying ones 
on the other side. This corresponds physically to an infinite sheet of 
supercurrent confined to the domain wall with a linearly rising gauge 
potential on one side and Meissner suppression on the other. While it 
was shown that these static solutions satisfied the field equations, their 
stability was left unresolved. In this paper, we analyse the normal modes 
of perturbations of the static solutions to demonstrate their perturbative 
stability. 

1 Introduction 

Topological and non-topological defects are interesting classes of solutions to 
study for a large range of physical systems. They are frequently manifest in 
cosmological models described by classical relativistic field equations and show 
up in quantised systems as non-perturbative effects. The more specific kind 
of defect known as a domain wall or kink can act as an interface, separating 
two regions described by the same physics but with different boundary condi- 
tions. Condensed matter physics uses domain walls to model phase transitions 
and large scale structures in a system. Additionally, domain wall solutions are 
used as a basis for branc- world models El EI El El Cj j where our universe is 
embedded in a higher dimensional space. This space is described throughout 
by the same physical model, but different asymptotic vacuum behaviour in the 
higher dimension induces a kink defect to which our 3 + 1 dimensional world is 
confined. 

In one such brane-world toy model, Rozowsky, Volkas and Wali I] have found 
interesting solutions consisting of a pair of concentric domain walls coupled to 
a pair of U(l) gauge fields. Their solutions look physically like infinite sheets 
of supercurrent confined to the wall, producing a linearly increasing gauge po- 
tential on one side and Meissner suppression on the other. The purpose of this 
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paper is to demonstrate the perturbative stability of this configuration. Our 
method is based on an analysis of the normal modes of perturbations of the 
static solutions, and we find that for a large range of parameters these modes 
are oscillatory, remaining bounded in time. 

In Section [3 we present the model with a slight generalisation and display 
the static solutions for light- and space-like gauge fields. Our stability analysis 
method is outlined in Section followed by a full investigation demonstrating 
the perturbative stability of the static configurations. We give a summary of 
our results in Section 0] 

2 The model 

We present a slight generalisation of the Rozowsky et al. model which specifies 
two scalar fields fa (i = 1, 2) each with independent local £7 (1) gauge symmetries 
with associated gauge fields A%. There is an additional discrete Z 2 symmetry 
interchanging fa <-» fa and <-> A%. This model is a toy model invented 
to study the clash-of-symmetries mechanism in its simplest non-trivial setting; 
see [HI El S3 ID] ■ A quartic scalar potential couples the scalar fields to each 
other and permits domain wall solutions asymptoting to different degenerate 
minima. The Lagrangian density is 

C = (D 1 * fa)* Dp fa + {D^fafD^fa - \F» v F lilv - \F£ v F 2llv - V(fa,fa), (1) 

with the scalar field potential given by 

V(fa,fa) = \i{falfa + fa 2 fa - v 2 ) 2 + X 2 fafafa 2 fa. (2) 

We work in the Ai ; 2 > parameter regime, where the degenerate global minima 
are manifestly given by 

\<f>l\ = v, 02=0 and fa — 0, \fa\ = v. 

The £7(1) gauge fields A^ are described in the usual way by F^ v and their 
appearance in the covariant derivative 

D» = d< l -iQ 1 At-iQ 2 A 2 l , 

where Qi is the charge operator associated with the £7(1), symmetry. Keeping 
the discrete 7L 2 symmetry, the £7(l)i <8 £7(1)2 charges of the scalar fields are 
fa ~ (e, e) and fa ~ (e, e), with e and e constants. 1 

The Euler-Lagrange equations of motion for the scalar and gauge fields are 

D^D^fa = -2Ai&(#& + fa, fa - v 2 ) - \ 2 fafajfa, 
d v F^ = 2Im(efaD>*fa + e<f>* fa ) , 

where i — 1 and j — 2, or i — 2 and j = 1 (this notation is to be un- 
derstood in subsequent equations). Following [Tj, we look for static solutions 
that depend only on z and utilise a polar decomposition for the scalar fields; 
fa{z) = Ri(z)e 1&i ( z \ The scalar potential (J2J allows one to construct domain 

lr The Rozowsky et al. model took e = 0, so this is a slight generalisation, first introduced 
in [Q] 
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wall solutions by requiring fa to asymptote to different degenerate minima of 
V in an appropriate limit. Since our fields vary only with z, this will be the 
dimension perpendicular to the wall and the limit will be as \z\ tends to infinity. 
This gives us the boundary conditions 

\fa\->0, \fa\^v as z -> -oo, ^ 



v, Ifal — > as z — > oo, 



or vice- versa. 

Since we are working with gauge fields we have the freedom to choose two 
gauges, one for each Af; the Lorentz gauge <9 M Af = turns out to be the most 
suitable choice for both. The algebra simplifies if instead of A? one considers 
the linear combination Af = eAf + eAj. Working with these choices, the field 
equations of motion reduce to 

R'l = -Ri(Af - Af - Af) + 2X 1 R. l (R 2 + R 2 - v 2 ) + \ 2 R t R 2 , (5) 
A (t,*,yf = 2(e 2 + ^ )R 2 A {t, X ,y) + 4e g jR 2_ 4 (t,x, S /) ; (6) 

Af = 0, (7) 
e< = -Af, (8) 

where prime denotes differentiation with respect to z. We immediately see that 
the A\ and hence the A\ are pure gauge and do not contribute to the physics; 
neither do the 6^. 

To further simplify the problem, we note that each gauge component Af ,x,y ^ 
exhibits the same dynamics in © and appears quadratically in JSJ. Thus, the 
qualitative physical behaviour depends on whether the gauge field configuration 
is space-like, time-like or light-like. Expressing this behaviour as the single field 
Ai we have 

R'l = kR t Af + 2\ 1 R l {R 2 l + R) - v 2 ) + \ 2 R,R 2 , (9) 
A'l = 2(e 2 + e 2 )R 2 Ai + 4eiR 2 Aj, (10) 

where k = +1, 0, —1 for space, light- and time-like gauge fields respectively. 
One can see that equation © is consistent with the boundary conditions J3J so 
long as k > 0. For k = — 1 the asymptotic behaviour of Rt is oscillatory and so 
we discard this time-like scenario. 

Considering equation (|1(J|) on the side of the wall where Ri — > v and Rj — > 0, 
we see that A" — > 2(e 2 + e 2 )v 2 Ai- Since physical solutions must be bounded, 
we conclude that Ai is exponentially suppressed. On the other side of the wall, 
Ri — > 0, Rj — > v and, using the previous result, Aj — > 0. Thus A" — > and 
this gauge field takes on a linear form. We note that all measurable quantities 
associated with Ai arise through derivatives and so this unbounded solution is 
still physical. 

The set of equations I© and l|10[l cannot, in general, be solved analytically 
and we use the relaxation-on-a-mesh technique to obtain numerical solutions. 

In the light-like case, k = 0, and Ri can be solved for independently of Ai. 
Typical solutions are shown in the top two plots in Figure The scalar fields 
assume a typical domain wall configuration asymptoting to distinct minima of 
the potential. As the boundary conditions for the Ri are symmetric under z 
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Figure 1: Static solutions for the two scalar and two gauge fields in the U(l) Cg> Z7(l) 
model, plotted against z. The top two plots are for the light-like case, the bottom 
two for the space-like case. All plots have e = 1, e = |, Ai = 1, A2 = 2 and v = 1. 
The plots on the left have symmetric boundary conditions for the gauge fields, those 
on the right asymmetric. In the light-like case, the scalar fields do not feel the gauge 
fields and thus do not depend on the choice of gauge field boundary conditions. This 
is unlike the space-like case where the scalar fields centre on the gauge fields to restore 
the reflection symmetry. 

reflection, the solutions for these scalar fields are just reflections of each other. In 
the left plot, the boundary conditions for the two gauge fields are also reflection 
symmetric. In the right plot, a different boundary condition is used for Ai- As 
the scalar fields do not feel the presence of the gauge fields they have exactly 
the same solutions in both cases. 

For the space-like case, k = 1, and the scalar and gauge fields are fully 
coupled. Solutions are shown in the bottom two plots in Figure □ with all pa- 
rameters, except k, mimicking the top two plots. Although they look similar, 
the light- and space- like plots on the left are slightly different. A more signifi- 
cant difference between these two scenarios is evident in the right plots where 
the boundary conditions for the two gauge fields are different. In the space- 
like case the scalar fields are influenced by the gauge fields and the favourable 
configuration is that with exact reflection symmetry. The boundary conditions 
serve to simply shift the centre of the domain wall and the right plot on the 
bottom is an exact translation of the left plot. Our result is contrary to the 
claim in pp that asymmetric boundary conditions in the space-like case are not 
equivalent to spatial translations of the domain wall centre. 

Disregarding the technical details, the qualitative features of this U (l)(g>t/(l) 
model are the reflection symmetric scalar fields in a domain wall configuration 
and the partially suppressed gauge fields. This suppression of Ai under their 
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respective Ri is physically similar to the Meissner effect and serves to semi- 
localise the gauge fields. We make the physical interpretation of an infinite 
sheet of supercurrent confined to the wall, producing a constant magnetic field 
in the region opposite the suppression. 

While we have shown the existence of static solutions to the model described 
by 0, we have not established their stability. In the next section we demon- 
strate that under small perturbations, the solutions to equations 10 and ijTU|) 
are stable. 



3 Stability 

The essence of static solutions is their time independence, but a physical pro- 
cesses requires the underlying fields to evolve in time. We must thus ensure 
that the static configurations found in the previous section are not destroyed 
by time-dependent perturbations. In this section we add to the static solu- 
tions perturbations expressed as normal modes and arrive at a set of equations 
characterising these normal eigenfunctions and associated eigenvalues. We then 
demonstrate that the eigenvalues are all positive and hence the perturbations 
are oscillatory. 

We begin by taking each static field, including all four gauge components 
A% , and adding a perturbation factored as an unspecified spatial part and a 
time dependent complex exponential. This exponential represents an arbitrary 
normal mode of the perturbation, characterised by an eigenvalue which is in 
general complex. We express this construction in the substitutions 

Ri ->i? 2 (z)+r l (z)e i ^*, 

At^At{z) + at{z)^\ (11) 

e 4 ^e 4 (z) + ^(z)e iwet . 

By the choice of an explicitly complex exponential, if uj is purely real then 
the perturbation will be oscillatory and hence remain bounded in time. On the 
other hand, if u> has an imaginary component, the exponential will blow up, 
signifying instability of the original static solution. 

We take the original field equations © , make the substitutions given by Ijllj) 
and simplify using the equations H, ||, @ and (JBJl for the static fields. We 
work to first order in r^, and 9i and consider only independent perturbations, 
which decouples the resulting set of equations to give 

(-flg - (Af - Af - Af) + 2Ai(3i?2 + R* - v 2 ) + A 2 i? 2 ) n = u^n, (12) 

(-5 2 2 + 2( e 2 + g 2 )i? 4 2 )af = cX, (13) 

(-%-2j±d^0 i =<Jt6 i . (14) 

Before we continue with these equations, we must first establish a general 
result. Given the equation 

f"(z) + V(z)f(z) + W(z)f(z) = 0, 

one can show that if W(z) < for all z, then there exist no non-trivial solutions 
for f(z) on the domain z £ R with f(z) — > as \z\ — > oo. To see this consider 
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large —z with / taking a vanishingly small positive value. For non-trivial so- 
lutions, / must increase as z increases 2 and so /' > 0. For solutions where / 
becomes vanishingly small for large z, we require /' < for some subsequent re- 
gion of the z-axis. This change in the sign of the first derivative requires /" < 
for some region, in particular we must have /" < when /' = 0, i.e. at the 
turning point. But at this point we have /" = —W(z)f and since / > and 
W(z) < for all z we have /" > 0. Thus the function is positive with a posi- 
tive gradient and can never turn back towards the z-axis. A similar argument 
holds when / is below the axis; it can never turn back up. Hence there are no 
non-trivial bounded solutions if W(z) < for all z. 

We now return to the issue of stability. Consider equation l|14(l with f(z) — 
0i(z) and W(z) = uj 2 ,. If oj 2 < then one would have W(z) < for all z and by 
the previous result the only solution for 9i would be the trivial one. Thus there 
are no negative eigenvalues for equation (|14J) with bounded eigenfunctions 9i. 
Note that the condition that eigenfunctions Qi be bounded does not preclude 
the analysis of the bounded nature of the perturbations. The perturbation to 
the static field is given in full by 0i(z)e luJl>t where by definition of a perturbation, 
9i(z) must be small and bounded. It is the nature of the temporal part e lu>at , 
hence the eigenvalues, that determines if the fields are stable. Since we have 
shown that w 2 > we have ojg real and thus oscillatory perturbations and hence 
a stable static field 0;. 

For the gauge fields, inspection of equation (|T3)l yields 

W(z) = u 2 a - 2(e 2 + e 2 )R 2 . 

For bounded af (z) we require W(z) > for some non-zero domain of z. This 
means that we need 

u 2 a > 2(e 2 + i 2 )R 2 

for some z. Since lo 2 is a constant it must be greater than or equal to the 
minimum of 2(e 2 + e 2 )R 2 , hence non-negative. Thus we have shown that the 
static gauge fields are stable under small time dependent perturbations. 

Following a similar argument for the scalar fields, equation l|12|) gives us the 
bound on the eigenvalues as 

w r > min(J7 (z)), 

where 

U(z) = -(Af - Af - Af) + 2Ai(3i? 2 + R 2 - v 2 ) + \ 2 R 2 . (15) 

It is not so clear as to the sign of this function. We analyse the light-like case first 
where the A^ terms are absent. In this case, as can be seen from equation 
with k = 0, the scalar field configuration and hence U(z) depend only on the 
parameters Ai, A2 and v. Since v can be absorbed into a rescaling of the we 
only have two parameters to consider. A typical plot of the function U (z) for 
the two permutations of i and j is shown on the left in Figure |2| 

We see that U(z) > for this choice of parameters. Figure shows the 
minimum of U(z) for a large range of values of Ai and A2. Since all minima 
are positive, it must be that w 2 > and hence the static scalar fields in the 
light-like case are stable, at least for this range of parameters. 

2 Since /(— 00) = there must be a region where / increases if it is to attain a finite positive 
value. 
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Figure 2: The function U(z) used to determine the eigenvalues of the Ti perturba- 
tion in the light- (left) and space- (right) like cases, plotted as a function of z. The 
parameters and corresponding field configurations are as in the reflection symmetric 
cases in Figure There are two plots in each graph corresponding to U(z) with i = 1, 
j = 2 and i = 2, j = 1. 
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Figure 3: The minimum of the function U(z) in the light-like case plotted against 
A2. The upper curves correspond to successively larger values of Ai, which runs from 
0.2 to 2 in steps of 0.2. 
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Figure 4: The minimum of the function U(z) in the space-like case. The plot on the 
left is against A2 with upper curves corresponding to larger values of Ai, which runs 
from 0.2 to 2 in steps of 0.2. The plot on the right is against the boundary condition 
for the gauge field with upper curves corresponding to larger values of v, which runs 
from 0.2 to 0.8 in steps of 0.1. 

In the space-like scenario, the results are similar. The plot in the right of 
Figure El shows equation (fT5|) with Af and A\ present. Figure 0] shows the 
minimum of U{z) for various values of the parameters and boundary conditions 
for the gauge fields. It is clear that the minima are all positive and so the static 
scalar field configuration is also stable in the space-like case. 

Since each field in the model permits static solutions which are independently 
stable, we conclude that the static configuration as a whole is a stable one. 
We have also verified this analysis with explicit numerical calculation of the 
eigenvalues. 

4 Conclusion 

Static solutions to a [7(1 ) ® U(l) gauged scalar model were recently found by 
Rozowsky, Volkas and Wali pQ. In this paper we have generalised the model 
slightly, presented the static solutions and demonstrated the stability of this 
field configuration. We achieved this by adding small time dependent perturba- 
tions, in the form of normal modes, to the static fields and obtaining eigenvalue 
equations. It was shown that these eigenvalues, corresponding to the normal 
modes, were positive for a large range of parameters in the model. Thus the 
perturbations were oscillatory and the static fields stable. 

Acknowledgments 

DPG was supported by the Puzey bequest to the University of Melbourne and 
RRV by the Australian Research Council. 

References 

[1] J.S. Rozowsky, R.R. Volkas, and K.C. Wali. Domain wall solutions with 
Abelian gauge fields. Phys. Lett. B, 580:249-256, 2004. 



8 



[2] N. Arkani-Hamed, S. Dimopoulos, and G.R. Dvali. The Hierarchy Problem 
and New Dimensions at a Millimeter. Phys. Lett. B, 429:263, 1998. hep- 
ph/9803315. 

[3] L. Randall and R. Sundrum. A Large Mass Hierarchy from a Small Extra 
Dimension. Phys. Rev. Lett, 83:3370-3373, 1999. hep-ph/9905221. 

[4] L. Randall and R. Sundrum. An Alternative to Compactification. Phys. 
Rev. Lett, 83:4690-4693, 1999. hcp-th/9906064. 

[5] K. Akama. Pregeometry. Lecture Notes in Physics, Berlin Springer Verlag, 
176:267-271, 1983. 

[6] V.A. Rubakov and Shaposhnikov M.E. Do we live inside a domain wall? 
Phys. Lett. B, 125:136-138, 1983. 

[7] I. Antoniadis. A possible new dimension at a few TeV. Phys. Lett. B, 
246:377-384, 1990. 

[8] A. Davidson, B.F. Toner, R.R. Volkas, and K.C. Wali. Clash of symmetries 
on the brane. Phys. Rev. D, 65(125013), 2002. hep-th/0202042. 

[9] E.M. Shin and R.R. Volkas. O(10) kinks: clash of symmetries on the brane 
and the gauge hierarchy problem. Phys. Rev. D, 69(045010), 2004. 

[10] A. Demaria and R.R. Volkas. Kink-induced symmetry breaking patterns in 
branc-world S77(3) 3 trinification models. Phys. Rev. D, 71(105011), 2005. 

[11] G. Dando, A. Davidson, D.P. George, R.R. Volkas, and K.C. Wali. Clash of 
symmetries in a Randall-Sundrum-likc spacetime. Phys. Rev. D. (in press, 
2005). 



9 



